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ABSTRACT 


Results  of  an  investigation  into  the  sensitivity  to  two  types 
of  inspection  error  of  link  sampling  procedures  described  by 
Harishchandra  and  Srivenkataramana  are  reported.  Relevant  compound 
distributions  are  derived.  Some  comparisons  with  results  obtained 
in  a  similar  investigation  for  standard  double  sampling  are  also 
given. 


Key  Words  and  Phrases:  Quality  control;  Compound  distributions; 
Binomial  distribution;  Hypergeometric  distribution;  Multivariate 
hypergeometric  distribution;  Two-stage  sampling. 


1 .  Introduction 


In  two  stage  acceptance  sampling,  it  is  necessary  to  inspect  further 
items  from  a  lot  when  the  evidence  provided  by  the  first  sample  is  incon¬ 
clusive.  Commonly  the  second  sample  is  twice  as  large  as  the  first,  though 
this  is  by  no  means  universal.  In  order  to  save  sampling  costs,  sometimes 
the  additional  evidence  is  provided  from  the  results  of  standard  ('first 
stage')  inspection  on  the  immediately  preceding  and  following  lots.  There 
is,  of  course,  an  implicit  assumption  that  these  neighboring  lots  are  of 
about  the  same  quality  as  the  lot  under  examination! 

In  this  paper  we  will  derive  some  formulas  for  the  properties  of  this 
kind  of  acceptance  sampling  procedure,  for  inspection  by  attributes  when 
inspection  is  not  perfect,  there  being  a  probability,  p,  of  (correctly) 
declaring  a  defective  item  to  be  defective,  and  a  probability,  p',  of 
(incorrectly)  declaring  a  nondefective  item  to  be  defective. 

2.  Link  Sampling 

A  procedure  described  by  Harishchandra  and  Srivenkataramana  (1982) 
(referred  to  as  HS  in  the  sequel)  is  as  follows: 

Routine  sampling  takes  random  samples  of  size  n  (without  replacement) 
from  lots  of  size  N.  Denoting  the  number  of  items  classified  as  defective 
(whether  correctly  or  not)  in  the  i-th  lot  by  Z.; 

if  <  a-|  the  lot  is  accepted 
if  Z^  >  a2  the  lot  is  rejected 

if  a1  <  Z1  i  a2,  the  quantity  Z^+Z^Z^  -  that  is,  the  total  number 
of  Items  found  defective  in  the  random  samples  of  size  n  from  the 
(i-l)-th,  i-th,  and  (1+l)-th  lots  -  is  calculated,  and 
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if  Z^+Z.+Z^  <.  &2  ,  the  lot  is  accepted 
lf  Zi-1+Zi+Zi+1  >  a2  ,  the  lot  is  rejected. 

The  numbers  a^.ag  and  a£  are  integers  chosen  to  give  some  desired  probabil Ities 
of  acceptance  for  specified  numbers  of  defectives  among  the  N  items  in  the 
i-th  lot.  Commonly,  but  not  necessarily,  aj  =  aJ,. 

3.  Distributions 

The  variables  Z..  j .Z^.Z^j  are  mutually  independent.  Generally 

Z.'-'(Bin(Y.,p)  *  Bin(n-Y.,p‘))  A  Hypg(n,D.,N)  (1) 

J  J  J  y  J 

j 

where means  "is  distributed  as"  ; 

*  denotes  "convolution"1 

A  is  a  "mixing"  or  "compounding"  symbol  indicating  the  Y  has  the 
Y 

distribution  following  the  symbol; 

D,  is  the  number  of  defective  items  in  the  j-th  lot; 

J 

Y.  is  the  number  of  (really)  defective  items  in  the  sample  of  size  n 

J 

from  the  j-th  lot; 

Bin(g,h)  denotes  the  binomial  distribution 

Pr[X*x]  =  (9)hx(l-h)9-x  (0  <  h  <  1;  x=0,l , . . . ,g) ; 

and  Hypg  (g,h,k)  denotes  the  hypergeometric  distribution 

Pr[X=x]  =  (jj)(g'jj)/(g)  (max(0,g-k+h)  <xi  min(g.h)) 

4.  Calculation  of  Acceptance  Probabilities 

From  (1)  we  can  obtain  the  explicit  expression 
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pru, .  i] .  i  k°J)(n'°J)  i  (j)(;:j)pup'2-uo-p)y-uo-p')"-p-z+" 
J  (JJ)  y  3  3  ua0  u  4  u 


(2) 


=  P(z|Dj) ,  say. 


(0  <  z  <  n) 

The  probability  of  acceptance  of  the  i-th  lot  is 

PrUi  £  +  P'"[(a1  <  Z.  <  a2)  (Zi_l+Zi+Zi+1  <  a£)] 


*  I  PUjlDj)  +  l  P(z1l01)  ,  ll  ,  /h-llOl-ll^wl^l)  (3) 

zi-l+21+l— a2"zi 


2f<a, 


al<zi— a2 


For  the  link  sampling  method  to  be  useful,  it  is  necessary  that  D^-j, 

D..  and  Di+1  do  not  differ  too  greatly.  In  the  case  of  binomial  sampling 
(corresponding  to  N  -*■  «,  the  case  considered  in  HS),  _if  the  lot  proportion 

defective  is  the  same  in  all  three  lots,  the  acceptance  probabilities  for 

link  sampling  are  the  same  as  they  would  be  for  regular  two-stage  sampling 

(with  the  same  values  of  8^2  and  a£)  with  the  second  sample  (of  size  2n) 

being  chosen,  when  needed,  from  the  i-th  lot.  A  similar  result  is  not  valid 
when  N  is  finite  (as  it  is  in  this  paper),  even  when  _ ^=D -=D (=D,  say) 
because  the  convolution 

Hypg  (n,D,N)  *  Hypg  (n,D,N)  (4) 

is  not  the  same  as  the  distribution 

Hypg.  (2n,D,N) 

(or,  indeed,  as  Hypg  (2n,2D,2N).) 

5.  Partial  Link  Sampling  Procedures 

It  may  sometimes  be  a  drawback,  in  the  link  sampling  procedure,  that 
It  Is  necessary  to  wait  for  the  results  of  inspection  of  the  (i+l)-th  lot 
before  reaching  a  decision  on  the  1-th  lot.  If  the  sample  from  the  1-th 
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lot  is  inconclusive.  One  way  around  this  difficulty  would  be  to  replace 

by  ( _1+Zi ) ,  which  can  be  calculated  immediately.  It  may  be 
felt  that  this  is  straining  the  assumption  of  (roughly)  constant  D  values 
unduly.  In  order  to  meet  this  difficulty  HS  propose  the  use  of  a  partial 
linking  sampling  procedure  in  which  a  second  sample  of  size  n  (not  2n)  is 
taken  from  the  i-th  lot  and  used  in  place  of  the  sample  from  the  (i+l)-th 
lot  is  reaching  a  decision.  This  means  that  (Z^ is  replaced  by 
(Zi  j +Zi+Z|),  where  ZI  denotes  the  number  of  items  found  to  be  defective 
(whether  correctly  or  not)  is  the  second  sample  from  the  i-th  lot.  The 
acceptance  probability  is 

PrCZx  aj]  +  Pr[(a1  <  Zi  <  a2)  n  (Z.^+Zj  <  a£)]  (5) 


Now,  of  course,  Z .  and  Z}  are  not  independent,  though  Z._j  is  independent 
of  (Z^.Z’. ).  The  joint  distribution  of  Z^  and  ZI,  obtained  from  equation 
(4)  of  Kotz  and  Johnson  (1983)  by  putting  nj  =  n2  =  n,  is 


Zi 


Bin(Y1,p)  *8in(n-Y1,p') 
Bin(Y2,p)  *Bin(n-Y2,p')jY1,Y2 


A.  Mult.  Hypg(D;n,n;N) 


(6) 


Here  Mult.  Hypg  (D;n,n;N)  is  a  multivariate  hypergeometric  distribution,  with 

(7) 

(0  yry2  i  n;  D-N+2n  <  yj+y2  i  °) 


» <v»2>]  - i-) 


(  "  ) 
'n,n,N-2n' 


Nl 


C(n)ir(N-2n)l 


The  expected  number  of  items  inspected  in  the  1-th  lot  is 
n{l+Pr[a^  <  Z^  <.  a2]}  ■  i[n+n{l+2  Pr[aj  <  Zj  <  «2]>] 
while  with  regular  two-stage  sampling  (with  second  sample  size  2n)  It  is. 
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n{l+2Pr[a.|  <  Z..  s  ag]) 

[HS  also  describe  another  modification  in  which  Z!  replaces  ^  rather 
than  Z^.  The  analysis  is  exactly  similar  to  that  set  out  above]. 

6.  Tables 

Table  1  gives  acceptance  probabilities  for  both  link  sampling  and  partial 
link  sampling  for  a1  *  1 ,  a^  =  a£  *  5  (the  values  used  in  HS);  and  N  =  100, 
200;  n  =  20,  50;  p  =  1,  0.9.  0.75;  p'  =  0,  0.1  and  four  sets  of  values  of 
proportions  defective 


WN 

D^N 

WN 

0.05 

0.05 

0.05 

0.1 

0.1 

0.1 

0.05 

0.1 

0.15 

0.15 

0.1 

0.15 

These  tables  were  computed  by  RNR.  The  multivariate  hypergeometric  distri¬ 
bution  in  (7)  can  be  computed  by  calculating  its  logarithm,  using  the  fact 
that  £n(k!)  =  £nr(k+l)  and  calling  a  log-gamma  subroutine  repeatedly.  However, 
this  can  result  in  a  severe  loss  of  significant  digits,  even  when  the  log-gamma 
function  is  evaluated  in  double  precision. 

A  better  alternative  (used  in  the  construction  of  Table  1)  is  to  re-express 
(7)  as  the  product  of  two  univariate  hypergeometric  probabilities: 
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The  equivalence  of  (7)  and  (8)  is  not  obvious,  but  it  is  easily  established 
algebraically  or  probabilistically.  The  advantage  of  computing  (8)  rather 
than  (7)  is  that  subroutines  for  the  univariate  hyperaeometric  distribution 
are  available  in  a  number  of  high-level  computing  languages. 

The  values  in  the  tables  exhibit  the  very  marked  influence  of  p'  on  the 
acceptance  probabilities. 

Table  2  contains  acceptance  probabilities  for  standard  double  sampling, 
with  the  same  parameter  values  (a^  =  1,  a2  =  a£  *  5)  as  in  Table  1.  The 
size  of  the  first  sample  is  n  ■  20  and  the  second  sample  (taken  from  the 
same  (i-th)  lot)  is  n"  =  2n  s  40.  The  probability  of  acceptance  is 


Pr[Z.  *  a13  +  PrC(a1  <  Z.  s  a')  n  (Z.+ZV  <  a2)] 


(9) 


where  ZV  denotes  the  number  of  item  judged  to  be  defective  in  the  second 
sample  (size  2n")  from  the  i-th  lot.  The  joint  distribution  of  Y^  and  YV, 
the  actual  numbers  defective  in  the  two  samples  is  the  multivariate  hyper- 
geometri c 


pr[(v*> n  vy-yjj  ■ 


[Given  Yi  and  YV,  Zi  and  ZV  are  distributed  independently  as  in  ( 2 .) 


(10) 

(10)’ 


The  figures  in  Table  2  should  be  compared  with  those  in  Table  1  for  n  *  20 
and  D^_.|/N=D,j/N=Di+i/N  (*  0.05,  0.1).  As  is  to  be  expected  the  partial 
link  sampling  acceptance  probabilities  fall  between  the  values  for  link  samp¬ 
ling  and  standard  double  sampling.  It  appears  that  they  are  closer  to  the  link 
sampling  values  than  to  the  standard  values.  The  differences  decrease  as  the 
lot  size  increases  (and  would  be  zero  for  infinite  lot  size). 
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Table  1:  Probabilities  of  Acceptance  for 


Link  Sampling  (LS)  and  Partial  Link  Sampling  (PLS) 


N  n 
100  20 


Acceptance  numbers  a^  *  1 ;  =  a£  *  5 

Probability  of  Acceptance 

On  1st  Sample 


Di-1/N 

D./N 

D.+i/N 

P 

P1 

LS 

PLS 

(LS  and  PLS) 

0.05 

0.05 

0.05 

1 

0 

.9572 

.9752 

.7395 

1 

0.1 

.2312 

.2302 

.1818 

0.9 

0 

.9733 

.9854 

.7802 

0.9 

0.1 

.2598 

.2591 

.2006 

0.75 

0 

.9887 

.9943 

.8379 

0.75 

0.1 

.3067 

.3065 

.2308 

0.1 

0.1 

0.1 

1 

0 

.5366 

.5464 

.3630 

1 

0.1 

.0808 

.0797 

.0735 

0.9 

0 

.6385 

.6535 

.4297 

0.9 

0.1 

.1043 

.1031 

.0922 

0.75 

0 

.7837 

.8012 

.5396 

0.75 

0.1 

.1504 

.1494 

.1269 

0.05 

0.1 

0.15 

1 

0 

.5348 

.6866 

.3630 

1 

0.1 

.0806 

.0880 

.0735 

0.9 

0 

.6373 

.7788 

.4297 

0.9 

0.1 

.1040 

.1144 

.0922 

0.75 

0 

.7835 

.8879 

.5396 

0.75 

0.1 

.1502 

.1652 

.1269 

0.15 

0.1 

0.15 

1 

0 

.4037 

.4547 

.3630 

1 

0.1 

.0746 

.0760 

.0735 

0.9 

0 

.4939 

.5585 

.4297 

0.9 

0.1 

.0948 

.0973 

.0922 

0.75 

0 

.6474 

.7202 

.5396 

0.75 

0.1 

.1343 

.1396 

.1269 

0.05 

0.05 

0.05 

1 

0 

.9503 

.9589 

.7372 

1 

0.1 

.2376 

.2371 

.1868 

0.9 

0 

.9680 

.9741 

.7763 

0.9 

0.1 

.2649 

.2645 

.2046 

0.75 

0 

.9857 

.9889 

.8327 

0.75 

0.1 

.3102 

.3101 

.2336 

0.1 

0.1 

0.1 

1 

0 

.5504 

.5550 

.3782 

1 

0.1 

.0872 

.0867 

.0789 

0.9 

0 

.6456 

.6522 

.4412 

0.9 

0.1 

.1102 

.1097 

.0972 

0.75 

0 

.7833 

.7910 

.5457 

0.75 

0.1 

.1553 

.1548 

.1308 
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Probability  of  Acceptance 

On  1st  Sample 


Di-1/N 

D./N 

WN 

P 

P‘ 

LS 

PLS 

(LS  and  PLS) 

0.05 

0.1 

0.15 

1 

0 

.5486 

.6782 

.3782 

1 

0.1 

.0870 

.0959 

.0789 

0.9 

0 

.6444 

.7637 

.4412 

0.9 

0.1 

.1099 

.1216 

.0972 

0.75 

0 

.7830 

.8717 

.5457 

0.75 

0.1 

.1551 

.1709 

.1308 

0.15 

0.1 

0.15 

1 

0 

.4225 

.4713 

.3782 

1 

0.1 

.0804 

.0822 

.0789 

0.9 

0 

.5081 

.5665 

.4412 

0.9 

0.1 

.1002 

.1032 

.0972 

0.75 

0 

.6538 

.7168 

.5457 

0.75 

0.1 

.1388 

.1447 

.1308 

0.05 

0.05 

0.05 

1 

0 

.2517 

.2041 

.1811 

1 

0.1 

.0018 

.0018 

.0018 

0.9 

0 

.3701 

.3543 

.2512 

0.9 

0.1 

.0029 

.0029 

.0029 

0.75 

0 

.5804 

.61 63 

.3786 

0.75 

0.1 

.0052 

.0052 

.0052 

0.1 

0.1 

0.1 

1 

0 

.0079 

.0078 

.0078 

1 

0.1 

.0001 

.0001 

.0001 

0.9 

0 

.0195 

.0191 

.0191 

0.9 

0.1 

.0002 

.0002 

.0002 

0.75 

0 

.0631 

.0594 

.0581 

0.75 

0.1 

.0007 

.0007 

.0007 

0.05 

0.1 

0.15 

1 

0 

.0079 

.0078 

.0078 

1 

0.1 

.0001 

.0001 

.0001 

0.9 

0 

.0195 

.0192 

.0191 

0.9 

0.1 

.0002 

.0002 

.0002 

0.75 

0 

.0629 

.0679 

.0581 

0.75 

0.1 

.0007 

.0007 

.0007 

0.15 

0.1 

0.15 

1 

0 

.0078 

.0078 

.0078 

1 

0.1 

.0001 

.0001 

.0001 

0.9 

0 

.0191 

.0191 

.01 91 

0.9 

0.1 

.0002 

.0002 

.0002 

0.75 

0 

.0582 

.0582 

.0581 

0.75 

0.1 

.0007 

.0007 

.0007 
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N 


200 


Probability  of  Acceptance 

On  1 s*  Samp! 


'i-l/N 

D./N 

Di+1/N 

P 

P' 

LS 

PLS 

(LS  and  PLS) 

0.05 

0.05 

0.05 

1 

0 

.3250 

.3198 

.2368 

1 

0.1 

.0028 

.0028 

.0028 

0.9 

0 

.4273 

.4290 

.2991 

. 

0.9 

0.1 

.0039 

.0039 

.0039 

0.75 

0 

.6051 

.6180 

.4108 

0.75 

0.1 

.0062 

.0062 

.0062 

0.1 

0.1 

0.1 

1 

0 

.0197 

.0195 

.0194 

1 

0.1 

.0002 

.0002 

.0002 

0.9 

0 

.0368 

.0361 

.0354 

0.9 

0.1 

.0004 

.0004 

.0004 

0.75 

0 

.0895 

.0872 

.0807 

0.75 

0.1 

.0010 

.0010 

.0010 

0.05 

0.1 

0.15 

1 

0 

.0197 

.0207 

.0194 

1 

0.1 

.0002 

.0002’ 

.0002 

0.9 

0 

.0367 

.0408 

.0354 

0.9 

0.1 

.0004 

.0004 

.0004 

0.75 

0 

.0892 

.1119 

.0807 

0.75 

0.1 

.0010 

.0010 

.0010 

0.15 

0.1 

0.15 

1 

0 

.0194 

.0194 

.0194 

1 

0.1 

.0002 

.0002 

.0002 

0.9 

0 

.0355 

.0355 

.0354 

0.9 

0.1 

.0004 

.0004 

.0004 

0.75 

0 

.0810 

.0819 

.0807 

0.75 

0.1 

.0010 

.0010 

.0010 
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Table 

N 

100 

0. 

0. 

0. 

0. 


_2: 

Acceptance 

Probabilities  With  Regular 

Double  Sampling 

In  all  cases 

n  =  20,  n* '=40;  &1  =  1#  ^ 

=a2 ' =5 

P 

* 

P 

Di/N  =  0.05 

Di/N  =  0.1 

1 

0 

1.0000 

0.5305 

1 

0.1 

0.2241 

0.0769 

9 

0 

1.0000 

0.6565 

9 

0.1 

0.2542 

0.0997 

75 

0 

1.0000 

0.8202 

75 

0.1 

0.3033 

0.1458 
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